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Abstract 

In view of the photon-number tomograms of two-mode light states, using the qubit- 
portrait method for studying the probabiHty distributions with infinite outputs, the sep- 
arability and entanglement detection of the states are studied. Examples of entangled 
Gaussian state and Schrodinger cat state are discussed. 
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1 Introduction 

The probability representation of quantum states with continuous variables was suggested in 
[1]. According to this representation, any quantum states is associated with a fair probability- 
distribution function called symplectic tomogram (or tomographic-probability distribution). In- 
formation contained in the tomographic-probability distribution is the same that is in the density 
operator. The optical tomograms |2l |3] were used earlier as a technical tool to reconstruct the 
Wigner function [4] which was identified with the quantum state. The experiments to measure 
quantum states [5] are aimed at obtaining the Wigner function by means of measuring the 
homodyne quadrature component. Though this procedure implies the measurement of optical 
tomograms, the latter ones were not identified with the quantum states in these experiments. 
The optical tomogram is considered as intermediate technical information that gives the oppor- 
tunity to "reconstruct" the state identified with the Wigner function. 

An understanding of the fact that tomograms (optical or symplectic ones) are themselves 
a primary notion of quantum states was suggested in [Ij (see also [SI El El E|)- The ana- 
logues treatment of quantum states of spin systems (qubits, qudits) as tomographic-probability 
distributions was proposed in jTOlITT]. Thus, in the probability representation of quantum me- 
chanics, which is completely equivalent to other representations like Schrodinger representation 
|12j . Feynmann path integral representation |13], Moyal representation [14j, etc., the quantum 
states are described by the standard probabiUty distributions (state tomograms). Spin tomo- 
grams were used to discuss the separability and entanglement phenomena in [T6l [T7] . The 
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qubit-portrait method for studying entanglement of multiqudit systems was suggested in [IB] . 
This method is similar to the integration method of symplectic tomograms applied in [19] (see 
also pO]), with both methods providing the analogs of two-qubit states. 

There exist descriptions of the photon states in terms of the photon-number tomograms 
[211 1221 123]. The photon- number tomogram is the probability distribution of a discrete variable 
n = 0,1,2,.... This distribution function also depends on an extra complex parameter and 
contains complete information on the quantum state. The photon-number tomograms for the 
Gaussian states were constructed in [23]. The aim of our work is to develop an analog of the 
qubit portrait method that was suggested for qudit states and apply this method to the photon- 
number tomogram of multimode photon systems. This means that we consider a linear map of 
discrete probability distributions with infinite probability vectors onto probability distributions 
with finite probability vectors. The map, which is applied to the photon-number tomograms, 
is not a positive map (see, e.g., p5]). Nevertheless, it provides a possibility to detect the 
entanglement of multimode photon states. Positive maps were used to detect entanglement in 
[2S1I271I2B] and nonpositive maps were used to detect entanglement, e.g., in [2^ [5U| ISTl IB^ IB5]. 

In this paper we focus on two- mode nonclassical states of light, which are of great interest 
since they can have some particular features [5^ [55] . Entanglement of quantum states 
turned out to be a very important tool for secure quantum communication and quantum cryp- 
tography. A geometrical interpretation of entangled two-qudit states was proposed in [36j- In 
spite of the fact that detection of entanglement has undergone rapid development in the last few 
decades, they do not cover all the possibilities. The main goal of the present article is to intro- 
duce an alternative method that can reveal the entanglement of two-mode light states. There 
have been already made successful attempts to solve this problem [2S1EZ], for instance, by using 
the Bell-CHSH inequality [2S1 [SH] within the framework of symplectic tomograms [H [T21 HP] . 
Here, we try to attack the problem with the help of the photon-number tomogram [2T| [22t [23] 
which provides additional information on quantum correlations. 

Since the entanglement of two qubit system has been sufficiently investigated, it looks feasible 
to detect the entanglement of the system involved by using a linear map of the photon-number 
tomogram with infinite outputs onto qubit tomogram (qubit portrait). The qubit-portrait 
method to study qudit states was introduced in [18] and developed in [15]. As far as there 
exist many ways to construct a qubit portrait, we are going to discuss some examples of them. 
The procedure above results in reducing the separability property of a two-mode light state 
to the Bell-CHSH inequality for two qubits. That inequality is fulfilled if the initial two-mode 
light state is separable. It naturally leads to a necessary condition of separability. Its violation 
indicates immediately that the state in question is entangled. 

The paper is organized as follows. 

In Sec. [2], we give a brief review of the photon-number tomogram and its generalization to the 
two-mode case. Here we also suggest a method of obtaining the joint probability distribution 
corresponding to a two-qubit tomogram. In Sec. [3], an example of Schrodinger cat state is 
studied in detail. In Sec. [31 we introduce the general approach to deal with the Gaussian 
states. A particular example of its application is given and another way of constructing the 
qubit portrait is demonstrated. In Sec. [5l conclusions and prospects are presented. 
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2 Photon-number tomogram and its qubit portrait 



The conventional photon-number tomogram for a one-mode hght state given by the density 
matrix p reads [2T| [22] [23] 

Tr {pb\a)\n){n\b{a)) = {n\b{a)pb^{a)\n), (1) 



win, a) 



where \n) is an eigenstate of the operator a^a, with a"'' and d being the photon creation and 
annihilation operators, respectively, and D{a) is a displacement operator, which depends on 
the complex number a as follows: 



b(a) = e"*^^- 



Thus, the photon-number tomogram is nothing else but the photon distribution function of 
the state {a)pD{a). In other words, it is the probability to find n photons in the state with 
the amplitude shifted by a complex number a. It is worth noting that once tomogram ([1]) is 
known, the state p can be reconstructed. (For Gaussian states the photon distributions were 
expressed in terms of Hermite polynomials of several variables in [HI W2\). 

As far as the two-mode light state is concerned, the formula ([1]) changes slightly 

w{ni,n2,ai,a2) = {nin2\b{ai,a2)pb\ai,a2)\nin2), (3) 

where \n1n2) = |^i)|^2) is the state with ni photons in the first mode and n2 photons in the 
second mode, and D{ai,a2) = Di{ai)D2{a2)- Here Di{ai) denotes the displacement operator 
([2]), where we replaced a aj, and d — Sj, with d| and dj being the photon creation 

and annihilation operators of the ith mode, i = 1,2. The operators Di{ai) and -02(^2) obviously 
commute. 

The photon-number tomogram ([1]) is rather similar to the spin tomogram flUi [TT] of a 
qudit (particle with spin j). The number of photons n plays the role of the spin projection m 
(m = —j, —j + 1, j), while the complex number a is analogues to the unitary matrix of the 
corresponding rotation group. In contrast to the spin tomogram, the photon-number tomogram 
has infinite outputs, so there are infinite ways to construct a qubit portrait [HI |15] of such 
a tomogram. Indeed, if one has the infinite probability-distribution vector with nonnegative 
components 



( ^(o,«) \ 

w{l, a) 
w{2, a) 

V ... J 



(4) 



where w{n,a) = 1, a new probability-distribution vector W2 with two components can be 

n=l 

constructed as follows 



( w(n, a) 

J2 w{n,a) 

\ ne(Z+\A) 



(5) 
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Figure 1: Partition of the set x Zj^ to subsets i = 1, ...,4, that enables one to construct 
the four- vector ([7]). Nevertheless, such a vector is not a probability distribution vector of two 
qubits because Bi cannot be written in the form ^ ^2 ■• where A^*"* and subsets of 

Z^. 



where A is an arbitrary subset of the set of nonnegative integers Z+ = {0, 1, 2, ...}. (We denote 
the set of nonnegative integers that do not belong to A as Zj\A] from this it follows that 

A u (z+\A) = z+, A n (z+\A) = 0). 

Let us now consider the photon-number tomogram of the two-mode light. 
Once the parameters a\ and 02 are given the function w(?7,i, 77-2, ai, ^2) can be treated as a 
table with an infinite number of both rows and columns 



( iy(0, 0, ai, 02) 1^(0,1,01,02) ^(0,2,01,02) ■■■\ 

1^(1,0,01,02) 1(7(1,1,01,02) w(l,2,Oi,02) ... 

1(7(2,0,01,02) 1(7(2,1,01,02) M7(2, 2, Oi, O2) ... 

V ••• ■■■ } 



(6) 



with the sum Yl, U7(ni, 71,2, oi, 02) being equal to unity. This distribution induces the four- 

n\ ,712=0 

dimensional probability distribution vector of the form 



W'a 



( E 1(7(^1, ^2, 01,02) \ 

(ni,n2)e-Bi 

E w(ni,n2, 01,02) 

(rai,n2)ei?2 

E 1(7(77,1,77,2,01,02) 

E 1(7(771,772,01,02) 

\ (ni,n2)e-B4 / 



(7) 
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where the sets Bj 



1,...,4) satisfy the conditions: Bi C x 



U B, 

i=l 



Zi X Z^, and 



BiH Bj = if z 7^ j. An example of separation of the set Z^ x Z^ that meets all the above 
requirements is illustrated in Fig. [H 

The probability vector obtained can be interpreted readily as the spin tomogram of a particle 
with spin j = 3/2, but it is hardly able to describe the two-qubit system. The matter is that the 
simply separable state of two qubits has the tomogram of a factorized form w{mi,m2, i^i, A^2) = 
w{mi, 1^ i)w{m2, '^2), where mi and m2 are the spin projections of the first and second spin-1/2 
particles to the axes A^i and 1^2, respectively. This feature underlies the Bell-CHSH inequality 
and is of a great importance. For the factorization above to be valid, it requires B^ to be the 
direct product of subsets of i = 1,...,4. Consequently, once purposing to construct the 
qubit portrait, one should rewrite formula ([7]) in the form 



/ E 722, 01,02) \ 

(?ii,n2) e A1XA2 

E w(ni,n2, 01,02) 

^ ^ ("1,12) e Aix{Z+\A2) 

E t(7(ni,n2, 01,02) 

(ni,7i2) G {Z+\Ai)xA2 

E w(ni,n2, 01,02) 

V (ni,n2) e {Z+\Ai)x{Z+\A2) 

where Ai and A2 are subsets of the set of nonnegative integers Z^ 
separable, i.e., 

w{ni, ^2, Oi, 02) = Wi{ni, ai)w2{n2, 02), 
then vector (IH]) can be written in factorized form 



(8) 



If the state is simply 



(9) 



/ E w{ni,ai) 

nieAi 

E w{ni,ai] 

\ nie{Z+\Ai) 



( E w{n2,a2) 

n2£A2 

E 1^(^2,02) 

\ n2&{Z+\A2) 



(10) 



which corresponds to the probability-distribution vector of the simply separable two-qubit state. 
As far as separable states are concerned, the four-vector ([H]) becomes merely the convex sum of 
the vectors (fTOj) . 

To draw an analogy with the probability-distribution vector of two qubits, we designate the 
components of the vector ([8]) as follows: 



:iii 



f w;(+,+, 01,02) \ 
w{+, -,01,02) 
ti?(-,+, 01,02) 
V w{-, -,01,02) J 

with all the components being nonnegative and their sum equal to unity. 

Comment. V The construction of the four-dimensional probability distribution vector from 
the matrix (jH]) can be presented in an elegant matrix forno. To start, we show how one can 



^Such a map for finite probability vectors was used in the work of M. A. Man'ko, V. I. Man'ko, and R. V. 
Mendes (unpublished) 
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obtain a probability-distribution vector III2 with two components with the help of a probability- 
distribution vector 



The general form of transformation Uf^ — > 11^2 reads 




(12) 



vf'2 




Pia + p2b + psc 
1 - {pia + p2b + psc) 



Pi P2 P3 

1 -pi I-P2 1 - Pa 




(13) 



where Pi > and = 1. Regarding 1^2 as a three-dimensional probability vector, one can 
write 





Pi 


P2 








I- Pi 


I-P2 




b 1 = MIZ?3, 


(14) 















where M is a stochastic matrix with no more than two nonzero rows. A generalization of this 
method to higher dimensions is obvious, and one has 







(15) 



V J 

To deal with the probability-distribution matrices, one can use their vector representation [18] 
or, equivalently, apply the technic above separately to the rows and columns. In the case of 
matrix we have 



/ 

w |_ w 









\ 



J 



(16) 



where the block 



provides a required qubit portrait. This completes the Com- 



\ w 

ment. A 

Since the four- vector flTTl) is merely a vector-function of parameters ai and 02, we can 
introduce new complex variables jSi and f52 and construct the stochastic matrix of the form 



M(ai,a2,/3i,/?2 



f w(+,+, 01,02) w{+,+,ai,(32) w{+,+,pi,a2) w{+, +, pi, P2) \ 

w(+,-, 01,02) w{+,-,ai,p2) w{+,-,pi,a2) w{+,-,pi,p2) 

w(-,+, 01,02) w(-,+,Oi,/32) w{-,+,pi,a2) w{-,+,pi,P2) 

V w(-, -,01,02) w{-,-,ai,p2) w(-, -,/3i,02) w{-,-,pi,p2) / 



(17) 
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The matrix f[T7|) is obtained just in the same way as its analog for a two-qubit system 
For this reason it exhibits the same properties. The most essential point is the following. If 
the initial state of two-mode light is separable then, the Bell-CHSH inequality is fulfilled for 
certain. That implies that the inequality 



Tr(M(ai,a2,/5i,/32) / 
holds true, where the matrix / is given by the formula 



< 2 



/ 1 
1 
1 

V -1 



1 \ 

1 

1 

-1 / 



(19) 



If there exist complex numbers ai, a2, Pi, P2 such that the inequality ( |T8l) is violated, 
the two-mode light state is entangled. Consequently, fulfilling the Bell-CHSH inequality is a 
necessary condition of separability of the two-mode light state. 



3 Detecting entanglement of the Schrodinger cat state 

To demonstrate the application of the technique described above, we consider the state that is en- 
tangled by construction. (One cannot help mentioning that the entanglement of the Schrodinger 
cat states was detected by means of symplectic tomograms in [20]). Here we investigate the 
Schrodinger cat state of the form B5] 

\^) = Ar(7i,72)(|7i,72) + 1-71,-72)), (20) 

where 171,72) = |7i)|72), states |7j) are eigenstates of the photon annihilation operators, i.e., 
O'illi) = 7i|7i), = 1,2 (coherent states of the first and the second modes, respectively). The 
factor iV(7i,72) can be found by employing the normalization requirement = 1. Indeed, 

since for arbitrary coherent states I7) and \S) the scalar product (7!^) can be written in the 
form 061 HZ] 



m 

it is of no difficulty to calculate 



expi --I7I 



(21) 



{iplij) = iV^(7i, 72) ((71,72171,72) +2Re((7i,72| -71,-72)) + (-71,-72! -71,-72)) 
= 2iV2(7i,72) (1 + e-2|^^l^-2l7.P) = 4iV2(7i,72)e-(l^^l'+l^^n cosh (|7i|' + I72I') ■ 

Therefore, 

exp{|(|7i|' + |72|')} 



^(71,72) 



2{cosh(|7i|' + |72|')} 



1/2- 



(22) 



(23) 
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In case of pure quantum state, i.e. p = \'^|J){^p\, the formula ([3D takes a more convenient form 



w^{ni,n2,ai,a2) = (nin2|-D(ai, 02)!^') 



(24) 



To get the exphcit expression of the photon-number tomogram in our particular case, we 
recall that I7) = D{'y)\0) and D{a)D{'y) = i)(a + 7)e("^*-"*^)/2_ xhis implies that the displace- 
ment operator transforms any coherent state into a coherent one. To be more precise, 

I)(«i,«2)|7i,72) = I)i(«i)|7i)^2(«2)|72) = e("^^^-"^^^+"^^5-"|7.)/2|c,, + + 72). (25) 

In much the same way, one can write 

^(«i,«2) I -71,-72) = A(«i)| -7i)^2(«2)| -72) = e(°^-^-"i^^+°^^^-"^^2)/2|^^ 

(26) 

The scalar product of the coherent state \6) and the Fock basis state \n) is 



5* 



{n\S) = {n\ie-^'^'/'Y.^\k)]=e 



-|5|V2_^ 



k=0 



k\ 



(27) 



Using the simplified formula and taking into account the results of Eqs. (1^ . (125]) . ([2SD, 
fl26l) . and fl27|) . we arrive at the following formula for the photon-number tomogram of the state 
under consideration 



Wsc{ni,n2,ai,a2 



-(|«lP + |a2|2 



X 



4ni!n2! cosh (|7i|2 + I72I2) 
("t7i+"272)(^^ + 7i)"H«2 + 72)"^ + e"^^^+"^^^(«i - 7i)"n«2 - 72)"^ 



f28) 



Now, when the photon-number tomogram is known, our goal is to make sure that the 
Bell-CHSH inequality can be violated, because the state involved is entangled. In Sec. [21 we 
emphasized that there exist many ways to construct the portrait corresponding to the two-qubit 
system. In our opinion, the simplest way is to split the tomogram components corresponding 
to the states: 

(i) the vacuum state, 

(ii) states with no photons in the first mode and with nonzero number of photons in the 
second one, 

(iii) states with nonzero number of photons in the first mode and with no photons in the 
second one, 

(iv) states with nonzero number of photons in both modes. 
Then the probability distribution vector ([8]) takes the form 



5c(ai,a2) 



w(0, 0, tti, ^2) 



E w(0,n2,ai,a2) 

n2=l 
00 

E w(ni, 0, ai, ^2) 

n\=l 
00 

E w(ni,n2,ai,a2) 

\ ni, 712=1 



( ^i'5C•(+, +, "1, ^2) \ 
w^sc(+,-,ai,a2) 
w^sc(-,+,ai,a2) 
V ti'sc(-,-,ai,a2) J 



(29) 
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O A1XA2 

^ AiX(Z+\A2) 
® (Z+\Ai)xA2 

O (Z,\Ai)x(Z,\A2) 



123456789 10 

Figure 2: Simplest example of dividing the set x into parts, that enables one to construct 
the qubit portrait of the photon-number tomogram. 



The pattern of such a partition is shown in Fig. O Using the explicit expression ( l28l) of the 
photon-number tomogram for our particular case, we get the components of of the four-vector 



wsci+,+,ai,a2 



2cosh(|7i|2+ |72|2) 



wsc{+,-,<yi,oi2) 



-(|aiP+|a2p 



cosh(^2Re(ai7i + ^272)) + cos(2Im(ai7i + ^272)) |, 

(30) 



2cosh(|7i|2 + I72I 



,|«2p+l72|- 



cosh(^2Re(Q;i7i) 



+ e 



|02p-|72|' 



cos(2Im(Q;i7i)) - cosh(^2Re(Q;i7i + ^272)) - cos(2Im(ai7i + 0:272)) k(31) 



-(|aiP+|a2|' 



wsc{-,+,ai,a2 
+ e 



2cosh(|7i|2 + |72|2 
l"iP-l7iP cos(2Im(a272)) - cosh(2Re(ai7i + 



el"iP+l7iP cosh(2Re(a272)^ 

«272)) - cos(2Im(ai7i + 0272)) k(32) 



wsci-, -, «!, 02) = 1 - wsc{+, +, Oil, 02) - wsci+, -, tti, «2) - Wsc{-, +, «!, 02). (33) 
Let us now consider the function of four complex variables (or eight real variables) 
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Figure 3: Violation of the Bell-CHSH inequality for the Schrodinger cat states of the form 
|7i! 72) + I — 7i5 —72)- Here the maximum value of the Bell number is denoted as /(71, 72). The 
qubit portrait is constructed using formula ( l29l) ("zero- nonzero" approach). 



Bsci(yi,a2,Pi,P2) 



Tr(M5C'(ai,a2,/3i,/?2 



where the matrix Msc{cii, ^2, Pi, P2) is given by the formula 



(34) 



Msc{o:i,a2,Pi,P2 



W^4 5c(ai,tt2) \^4Sc{ai,P2) 5c(/5i, "2) sc{Pi,p2) ) 



(35) 

and the matrix / is determined by Eq. ( fT9l) . 

To demonstrate that there exist such values of variables ai, a2, Pi, and P2 such that the 
value of function Bsc{(^i,C(2, Pi, P2) is greater than 2, we find the maximum of this function, 
which depends entirely on the form of the quantum state 



/(7i,72)= max Bsciai, ^2, Pi, P2)- 

ai,a2,Pi,P2SC 



(36) 



Numerical calculations show that, whatever values of complex variables 71 and 72 are taken, 
the value of function /(71, 72) is never less than 2. This fact is depicted in Fig. [31 Nevertheless, 
the extent to which 7(71,72) is greater than 2 differs from one state to the other. If |7i| ^ 1 
and I72I ^ 1, then the efficiency of the technique proposed essentially decreases. 

4 Detecting entanglement of the Gaussian states 



The next example to be investigated is a nonclassical state of light described by the Gaussian 
density matrix (see, e.g., [281 ESI [371 HH])- It is of some interest because it can characterize 



10 



a two-mode mixed squeezed state. A Gaussian density operator p of the two-mode light is 
described by the Wigner function W{pi,p2, Qi, (I2) of the form [H] 



where 



W{p^,p2,qi,q2) = -^=l=exp|-^(Q- (Q))M-i(Q 



(Q)) 



(37) 



Q 



V q2 / 



M 



( '^pipi 

Or-. 



'P2P1 



qrpi 



q2pi 



0", 



PlP2 



a. 



P2P2 



qiP2 



92 P2 



Pi 91 



P2gi 



cr. 



91 91 



a. 



9291 



Pi 92 



P292 



9192 



a. 



9292 



(3J 



Here averaging (A) imphes {A) = TipA, where A can be one of the operators pj = —i{aj — 
a])/v^, qj = {hj + at)/v^, j = 1,2 (variables are assumed to be dimensionless, with % = 1). 
The components of the real symmetric dispersion matrix M read 



c^P.P, = {ViVj) - {Pi){Pj), c^9.9, = im) - ^P.9, = ^gJP^ = ^{Viqj + ijPi) - {Pi){^j)- (39) 

Thus, all the information on the state is contained in the Wigner function ( 137|) depending 
on 14 real variables (four of them determine the vector (Q) and ten of them define the matrix 
M). 

To obtain the photon- number tomogram ([3]), it is worth noting that this is merely a photon 
distribution function of the state Z)(ai, a2)p-D^(ai, a2)- Let us now reveal what is the Wigner 
function of such a state. For this purpose, we calculate the new average values (see, e.g., [21] )• 
For instance. 



{Pi)aia2 = Tr yD{ai, a2)pD\ai, ^2)^1 

= Tr (pZ)I(ai)Z)J(«2) {-^(ai - a\)/V2} bi{ai)t)2{a2) 
= Tr (^p |— i(ai + tti — a| — aD/V^^^ 

= Tr (p [pi + Im ai}) = (pi) + V2 Im ai. (40) 

Here we used the relation t)"^ {a)aD{a) = d + al. In a similar way we find (^2)0102 — (^2) + 
\/2Ima2, (^1)0102 = (q'i)+-\/2 Re «!, and {q2)aia2 = (q'2)+\/2 Re a2- The remarkable fact is that 
the elements of the matrix M do not change through a transition from p to D{ai, a2)pD^{<yi, 0^2)- 
It follows directly from the formulas (|39ll . We will refer to the vector of new average values as 

(Q)aia2 ■ 

The photon distribution function of a generic A^-mode mixed Gaussian state of light has 
been obtained earlier [2HI1I1I12]. We use it to compute the photon- number tomogram ([3]). The 
result is 

, , exp {-(Q),,,,(2M + U)-\Q)a,a,} i^S2,ni.n2(y) 

wcini, n2, ai, ^2) = , = , , (41) 

/det fM + 114) 
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where I4 is the 4x4 identity matrix, and the four- dimensional matrix R and the four-vector y 
can be expressed with the help of the matrix 



U 



/ -i i 0\ 

1 -i i 

71 1 1 

V 10 1/ 



(42) 



as 



R = U"^ (I4 - 2M) (I4 + 2M)-^ U 
y = 2U*^ (I4 - 2M)-^ (Q)„, 



02 1 



(43) 
(44) 



and the term of the form H^^l^ ^.^ ^.^ (x) is the four-dimensional Hermite polynomial given by 
the formula HQ] 



rrm 

ki,k2,k3,k4 ' 



X) = (_l)fci+te+fc3+fc4 ( _xRx 



Qki+k2+ks+k4 



2 J dx'l'dx^2'dxl'dx^^'' 



-o E 



ijXiXj 



(45) 

Once one has the expression of the photon- number tomogram (HTl) . it seems possible to 
explore the entanglement of the Gaussian states just in the same way as in Sec. [31 By analogy 
with the previous example, we choose the state which is entangled for sure. For instance, the 
squeezed state given by its Wigner function of the form 



W{pi,p2, qi, 92) = 4exp |-2 (Spj - V35pip2 + Spj + ql - V3qiq2 + qf) } 



(46) 



is entangled, because it contains the products pip2 and qiq2 in the exponent. This state corre- 
sponds to the matrix M and the vector (Q) of the form 



M 



(Q) 



(47) 



3 V35/2 
^35/2 3 

1 V3/2 

V3/2 1 

which meets the generalized uncertainty relation det M > 1/4^. The substitution of M for such 
a matrix in Eqs. (H3|) and (jH]) yields 



V 



/ 






V / 



R 



1 
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3^35-7^3 



3^35-7^3 

-3^35-7^3 
-3v^- 7^3 



-3^-7v^\ 

-3^35-773 

3^35-7^3 

3^35-7^3 / 



(4J 



and 
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y = 




/| Re aa) 
/| Re ai) 
I Re 02) - 
I Re ai) - 



- \/3 Im a2 ^ 

- V3 Im 
-\/3 Im a2 

Im «! y 



(49) 



Now the photon- number tomogram fHT]) can be calculated for any given parameters rii, 122, 
«!, and 02- It is just the right time for trying to detect the entanglement of the state (H6l) . 
To start, we apply the approach developed above, which means constructing the probability- 
distribution vector (l29l) . and then the stochastic matrix (!35|) . the trace of the product of which 
with the matrix (fT9l) leads to the Bell-like number. The obstacle arises during the first step, 
because the probability-distribution four-vector (l29i) is hardly computable analytically. The 
numerical calculation has some particular features concerned with four-dimensional Hermite 
polynomials. The matter is that it is rather difficult to calculate H^\^ „^ (y) when rii ^ 1 
and ?T,2 ^ 1. Therefore, we deal with numbers rti and n2 confined by the ranges < rii < 30, 

< ?T,2 < 30. To neglect the values of Wc{ni^n2^ ai, 0:2) providing rii > 30 or n2 > 30 we impose 
the limitations |Reaj| < 2, |Imai| < 2, i = 1,2 (and consequently |Re/3i| < 2, |Im/5j| < 2, 

1 = 1,2). If this is the case, the total contribution of probabilities w(ni, n2, ai, ^2) with high 
photon numbers (ni > 30, n2 > 30) to the distribution vector (l29l) does not exceed 10~^ and 
serves as an error measure. The calculation shows that the Bell-like number (1341) is less than 2 
providing the variables ai, 02, Pi, and P2 are not beyond the region specified above. 

The problem occurs not due to restrictions introduced but because of the unlucky choice 
of the qubit-portrait construction. Indeed, during the maximization process the variables ai, 
a2, Pi, and P2 attempt to distribute in such a way that the quantity wg'('^i, '"'2, «i, ^2) is far 
from zero when rii ^ 722 ~ 0, the quantity WG{ni,n2,ai, (32) is far from zero when rii ^ 
and n2 ^ 1, the quantity Wg('T-i, "^25 A, 0^2) is far from zero when rii 3> 1 and n2 ~ 0, and 
the quantity WG{ni,n2, Pi, P2) is far from zero when rii ^ 1 and n2 ^ 1. In Fig. [2|, one can 



see that the narrowness of the domains B2 and B3 prevents the quantities w{0,n2,ai, P2) 



and £ w{ni,0, Pi,a2) from increasing. In our case they cannot be greater then 2/3. (See 



analogues problem discussed in [50] ) 

To put the situation right we suggest the other separation {Ai x A2, Ai x (Z^\A2), {Z+\Ai) x 
A2, {Z+\Ai) X (Z_|_\A2)} of the set Z+ x which is needed to construct the four- dimensional 
probability-distribution vector. The main idea of the partitioning proposed is patterned in Fig. 
m It implies distinguishing even and odd numbers of photons in individual modes. Such a 
choice is stimulated by the peculiarities of the photon-distribution function of squeezed states 
concerned with its oscillations [52t [5T] . Thus, the explicit expression of the four- vector M4 g 
reads 



00 



n2=l 



00 



ni=l 



13 




Figure 4: Alternative partition of the set x Zj^. The sets Ax and A2 are nonnegative even 
integers, while the sets Zj\A\ and Zj\A2 are nonnegative odd integers. 



00 00 



4 Glttl, a2j 



E E w(2ni, 2^2, 01,02) 

ni=0 rt2=0 



00 00 



E E w(2ni,2n2 + !,«!, 02) 

ni=0 n2=0 
00 00 

E E w(2ni + 1,2^2, 01,02) 

ni=0 n2=0 
00 00 

. E E w(2ni + 1,2^2 + 1,01,02) 

\ ni=0n2=0 



wg(+, -,01,02) 
wg(-,+, 01,02) 
V «^g(-, -,01,02) / 



(50) 



The matrix Mg(oi, 02, /3i, (i-i) appears naturally in the same way as the matrix (!35]) . To demon- 
strate that the inequality [Tr (Mg(oi, 02, /3i, /32)/)| < 2 can be violated, it is enough to make 
the substitution oi = — 0.12z, 02 = 0.04i, /3i = 0.22z, (^2 = — 0.32z. Then the matrix Mq takes 
the form 



Mai-O.Ui, 0.042, 0.22i, -0.32i) 



f 0.6199 0.5907 0.6083 0.4678 

0.0222 0.0515 0.0291 0.1696 

0.0241 0.0395 0.0357 0.1624 

V 0.3335 0.3181 0.3266 0.2000 / 



(51) 



and |Tr(MG(-0.12z,0.04i,0.22i, -0.32i)/)| ^ 2.26 > 2. The violation of the Bell-CHSH in- 
equality for this particular set of variables oi, 02, Pi, and P2 indicates immediately that the 
Gaussian state fH6l) is entangled. 



While dealing with Gaussian states, one should notice that the presented method of to- 
mogram calculation works well for both pure and mixed states. This is accompanied by the 
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1 .7 "I ' ' ' ' ' ' ' ' — 

0.5 1 1.5 2 2.5 3 3.5 A k 



Figure 5: Violation of the Bell-CHSH inequality for the family of states fl5^ . The maximum of 
the Bell-like number is denoted by f{k,l). The purity of states does not change while moving 
along the separate curve: pure states with / = (a), mixed states with / = 0.01 (b), / = 0.04 
(c), and / = 0.07 (d). 



complexity of calculations including multidimensional Hermite polynomial contrast to the rela- 
tively easy ones for pure states. To demonstrate the advantages of the technique proposed, we 
consider the family of generally mixed states given by matrix M and the vector Q of the form 



M{k,l) 



I 



\ 












k 











k 








Jk^ 








\ 



A;2 - 
k 



J 



(Q) 






V J 



The purity of this state is determined by the parameter / only, because det M 
Thus the state is pure when / = and mixed for all / > 0. Denoting 



(52) 



;i+4/)/i6. 



f{k,l)= max Tr(S(ai,a2, A,/?2)/) (53) 
oi,02,Pi,P2e(-' 

we succeeded in plotting the dependence of f{k, I) on k for the states with different purities (see 
Fig. [5]). Though the violation of the Bell-CHSH inequality is clearly observed for pure states, it 
is rather difficult to detect entanglement of mixed states. This fact is similar to the two-qubit 
system behaviour when the increase in mixing causes the loss of entanglement. 

We have been just convinced that the technique of constructing the qubit portrait proposed 
in this section can be efficient enough for detecting entanglement of Gaussian states. To check 
if it works properly for other kinds of states, we apply it to the Schrodinger cat state discussed 
in the previous section. Though the photon-number tomogram of such a state can be obtained 
by means of addition of probability distribution functions of Gaussian states, we are going to 
exploit the explicit formula fl28l) . 

Then the even-odd photon number form of the probability distribution vector fl50|) reads 
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= ^ |e-'^<"^^^+"^^^) cosh (|«i + 7iP) cosh (|«2 + 72p) 

^g2Re(at7i+a^72) ^^gl^ _ ^^|2) cosh (1^2 " 721') 

+2cos(2ImK7i + a^72)) 
X cosh (loiip - |7ip) cosh (|Q;2p - 172^) cos(2Im(Q;i7i)) cos(2Im(Q;272)) 
-sinh - |7ip) sinh (|a2p - I72P) sin(2Im(a^7i)) sin(2Im(Q;272)) 

+2sin(2Im(at7i + a^72)) 
X cosh (|q;i|^ - |7i|^) sinh (|q;2|' - I72I') cos(2Im(ci;*7i)) sm(^2Im(Q;272)) 

+ sinh(|Q;i|2 - |7i|2)cosh(|Q;2|' - I72P) sin(2Im(Q;*7i)) cos(2Im(Q;;72)) |, 
w{+, -, a„ a,) = ^XC^ltl) ^ |e-'^<'^^^^+"^^^) cosh {\a, + 71^) sinh {{a^ + 72P) 

^g2Re(a^7i+a572) ^osh - 7i ^ sinh {\a2 - 72!') 

+2cos(^2Im(a^7i + 0272)) 
X cosh (Ictip - |7i|2) sinh (|Q;2p - 172^) cos(2Im(Q;;[7i)) cos(2Im(Q;^72)) 

-sinh(|aip - |7ip) cosh(|«2p - I72P) sin(2Im(a^7i)) sin(2Im(a272)) 
+2sin(2Im(a*7i + a*72)) 
X coshdccip - |7ip)cosh(|Q;2p - 172^) cos(2Im(Q;^7i)) sin(2Im(Q;;72)) 

+ sinh(|Q;ip - |7ip) sinh(|Q;2p - 172!') sin(2Im(Q;i7i)) cos(2Im(Q;272)) 
w{-,+,a„a,) = :XC^ltl) ^ {e-'^<"^^^+'^^^^) sinh (|ai + 7iP) cosh (1^2 + 72I') 

^g2Re(at7i+a^72) gmh - 7i|2) COsh (|«2 - 72^) 

+2 cos(2Im(a*7i + a*72)) 
X sinh (Ictip - |7ip) cosh (|Q;2p - 1721') cos(2Im(Q;;f7i)) cos(2Im(Q;;72)) 

-cosh(|aip - |7ip)sinh(|a2|' - I72I') sin(2Im(ai7i)) sin(2Im(a272)) 
+2sin(2Im(a*7i + a*72)) 
X cosh (Iccip - |7ip) cosh (|Q;2p - 172!') sin(2Im(Q;^7i)) cos(2Im(Q;;72)) 



(54) 



(55) 
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Figure 6: Violation of the Bell-CHSH inequality for the Schrodinger cat states of the form 
|7i! 72) + I — 7i5 —72)- Here the maximum value of the Bell-like number is denoted by /(71, 72). 
The qubit portrait is constructed using formula ( l50l) ("even-odd" approach). 



+ sinh (|aip — |7in sinh (|a2p — 172^) cos(2Im(a]'7i)) sin(^2Im(a272) 



(56) 



w{-,-,ai,a2) = 4cofh(|7ipl|72n ^ 



-2Re(a*7i+«*72) gj^^ (|ai + 71^) siuh (1^2 + 72 | 



_^g2Re(a*7i+a272) gj^^ (|ai - 71^) siuh {\a2 - 72!') 



X 



+2 cos(2Im(a^7i + 0372 
2 



sinh (laip — |7iP) sinh (|a2p — |72l^) cos 



— cosh (|aip — I71P) cosh (|a2p — 172^) sin 



X 



(^2Im(a*7i)) cos(2Im(a272 
2 _ u,j2^ c5ij^|^2Im(a*7i)^ sin(^2Im(a272) 

i72: 



+2 sin(^2Im(a;*7i + 0:272 
sinh (laip — |7in cosh (|a2p — 172^) cos(2lm{al^i)^ sin(2lm{al 



-Fcosh(|ai|2 - |7i|2)sinh(|a2p - 172^) sin(2Im(a^7i)) cos(2Im(a^72)) |. (57) 

By a method analogues to that used in Sec. [3l we calculate the maximum of the Bell-like 
number Bsc{c(i,C(2, Pi, P2) = \Tr{Msc{o:i,a2, Pi, ^2)1)1 over the complex numbers ai, 02, Pi, 
and P2, which is a function depending entirely on the form of the quantum state 



/(7i,72)= max Bsc{ai, a2. Pi, p2)- 
"i,"2,pi,P2ec 



(5J 



17 



The numerical calculation of this function for real numbers 71 and 72 is shown in Fig. 
[6l Comparison of Fig. [6] with Fig. [3] apparently shows that the "even-odd" approach of 
constructing the qubit portrait achieves better results than the " zero- nonzero" method. 

In Fig. [6], we can see that the Schrodinger cat state (l20l) is always entangled except the cases 
7i = or 72 = 0, where it becomes separable |0)(|72) + | — 72)) or (I71) + | — 7i))|0). While 
moving to large |7i| and |7i|, the function /(7i,72) increases. This increase can be ascribed to 
the fact that if |7i| 3> 1, the states |7j) and | —72) become quasiorthogonal due to | (7i| — 7i)| -C 1. 
If this is the case, the state (12(1 resembles the usual two-qubit cat state -^(lOO) + |11)), which 

is maximally entangled and leads to the greatest violation of the Bell inequality B = 2^/2 (also 
known as the Cirelson bound 153J). As far as our case of two- mode light state is concerned, 
the function /(7i,72) equals approximately 2.7 when 71 = 72 = 10, equals approximately 2.78 
when 7i = 72 = 50, and seems to tend to 2^/2 when |7i| and I72I go to infinity. 

5 Conclusions 

The photon-number tomogram is similar to the qudit tomogram (i.e., spin tomogram of the 
system with the total spin j), because of the countability of the outputs. The conventional 
approach of constructing the qubit portrait of qudit states can be applied to the pnoton number 
tomogram with some modifications. Being as powerful as its qudit counterpart, such a method 
suits for detecting entanglement of two-mode light states. To be more exact, the violation of 
the Bell-CHSH inequality indicates immediately that the state is entangled. Since there are 
many ways of reducing the photon-number tomogram of the two-mode light state to the spin 
tomogram of the two-qubit system, the effectiveness of the entanglement detection depends on 
the method utilized. To illustrate this fact we compared two of them in this paper. 
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